Interferometry with Synthetic Gauge Fields 



Brandon M. Anderson, Jacob M. Taylor, and Victor M. Galitski 
Condensed Matter Theory Center and Joint Quantum Institute, 
Department of Physics, University of Maryland, College Park, MD 20742-4111 

We propose a compact atom interferometry scheme for measuring weak, time-dependent accelera- 
tions. Our proposal uses an ensemble of dilute trapped bosons with two internal states that couple to 
a synthetic gauge field with opposite charges. The trapped gauge field couples spin to momentum to 
allow time dependent accelerations to be continuously imparted on the internal states. We generalize 
this system to reduce noise and estimate the sensitivity of such a system to be S ~ 10~ 7 



In recent years atom interferometry has emerged as a 
powerful tool for precision gravimetry and accelerome- 
try. [TH3] Experiments such as [3HS] are the most ac- 
curate measurements to date of surface gravity, cer- 
tain fundamental constants [TfjHT2"] and also provide 
probes of General Relativity and the inverse square 
law [H [13l [14] . Furthermore, accelerometers have wide 
application in more practical settings such as inertial nav- 
igation, vibration detection, and gravitational anomalies 
such as oil fields [Jjj] . Current experiments use short Ra- 
man pulses to manipulate spin states followed by periods 
of free evolution, corresponding to free flight, to accu- 
mulate sensitivity to external fields. During free flight, a 
sensitivity to external fields is imparted on the internal 
spin states in the form of a path-dependent phase. This 
phase can then be measured through a final Raman pulse 
and spin-dependent fluorescence techniques. 

At the same time that interferometry has emerged as 
a tool, interest in synthetic gauge fields has also arisen, 
mostly in the context of the quantum Hall Effect [TBI HZ] 
and cold atom spintronics [TBH2T] . These systems use 
optical coupling of internal spin states, simultaneous with 
momentum exchange with Raman laser beams, to induce 
an effective vector potential. Depending on the optical 
configuration these setups can simulate systems such as 
spin-orbit coupling 22,23, monopoles [UJ , or a constant 
magnetic field [Ml 125] . 

The optical coupling to the internal degrees of freedom 
provides a continuous coupling of momentum and spin. 
This is in contrast to standard interferometry schemes 
where spin and momentum coupling is generated only 
through a set of discrete Raman 7r/2 and ir pulses. In 
this paper we propose a new type of interferometer that 
uses the spin-momentum coupling to measure AC signals. 
We use the continuous spin-momentum compuling of the 
gauge field to produce an interferometer sensitive to high 
frequency time dependent (or AC) fields. This is in con- 
trast to current systems that at best are sensitive to con- 
stant (DC) or low frequency < 1Hz signals [13 [57]. We 
specifically propose using a trapped system of cold bosons 
under the influence of an optically induced gauge field 
to measure weak AC gravity signals. We discuss some 
potential implementations and we estimate that such a 




FIG. 1. A potential implementation of our interferometer 
based upon Ref. [IT] . The Raman beams £ii,2 couple a three 
level atom by two parallel Gaussian profiles with peaks that 
are spatially offset. Two of the dressed states become degen- 
erate in the large detuning limit, A — >■ oo, and couple to a 
"synthetic gauge field" with opposite charges. 

system will have a sensitivity of S ~ 7 ^^=- . We note 
that since our system is trapped it can be implemented 
on an atom chip 28-30 . 

As a toy model, we consider a single particle with an 
internal degree of freedom (pseudo-spin) in a harmonic 
trap with spin-orbit coupling and an external force: 

H = [ P'^ (f ~ )]2 + - mg(t) • f (1) 

where p and f are the position and momentum opera- 
tors respectively, u>o is the trapping frequency, cr = ±1 
labels the pseudo-spin of the particle, m is the mass of 
the particle, g(t) is the time dependent external force and 
A(f ) is the spin-orbit coupling field, or vector potential. 
We confine the particle to a two dimensional plane and 
chose vector potential to have the form of a magnetic 
field A(f) = mui c xy, where w c is the characteristic fre- 
quency scale of the spin-orbit coupling and e z is the unit 
vector perpendicular to the plane of confinement. This 
toy model captures the ideal behavior of systems such as 
Ref. Q7l [25]. 

Without spin-orbit coupling, (cj c — > 0), the path of the 
particle will depend on the force g(t) and will be indepen- 
dent of spin. With spin-orbit fields the path will depend 
on the spin of the particle as well as the force. In an atom 
interferometer, path differences can be mapped to an in- 
terference signal by creating an initial superposition of 
two spin states and sending them on spin-dependent tra- 
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jectories. The phase picked up on a semi-classical path 
can be found by taking the variation of the action S in 
the path integral formulation of quantum mechanics. 

For the system in Eq. [T] we find the first order phase 
due to g is 



JS c i _ /dtr(t)-g(t) 



(2) 



We become sensitive to this phase by creating an ini- 
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, |f)|0i) +e' 5 - 1 |-l)|(?!>_i)). In general the spin 
states |±1) are entangled to the orbital states \4>±i). 
In order to measure the phase we desire a pure spin 
measurement. We chose the paths to have complete or- 
bital overlap, (0i|</>_i) = 1, at the time of phase mea- 
surement. This places the system in the pseudo-spin 
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T)) which allows us to measure 
AS = S-i — Si though a single operator measurement 
such as S y . 

The physics described in our toy model is an idealized 
version of the proposal given in Ref. [17] , although other 
setups such as the experiment by Y.J. Lin et. al [25] have 
similar physics. This setup uses a cold atom in A-scheme 
as can be seen in the inset of Fig. [I] The Raman beams 
used have a Gaussian profile with offset centers which 
give spatial dependence to the dressed states. This spa- 
tial dependence induces dynamics that is identical to a 
charged particle in a magnetic field. In the large detun- 
ing limit, A — » oo, one of the bright states becomes de- 
generate with the dark state, however the "charge" that 
these two states see has opposite sign. We note that the 
synthetic field in this setup is non-uniform. This leads 
to further technical difficulties but does not change the 
underlying physics of our proposal. 

We now detail the specific solution for our spin-orbit 
coupled system described by Eq. [T] We start by solving 
the Heisenberg equations of motion for the system. These 
correspond to the classical, spin-dependent, Hamilton 
equations of motion since the system is quadratic in r 
and p. Not surprisingly, these solutions correspond to 
a combination of cyclotron orbits and orbits around the 
trap center. The direction of the orbits is set by the 
pseudo-spin, a. 

The solutions are characterized by the two frequencies 
uj± = uj±uj c /2, with Cj 2 = ujq + (oj c /2) 2 . We will consider 
initial conditions given by r(0) = ro and r(t) = 0, as 
will be discussed later. Fig. |3(b)| shows these paths in 
the absence of a driving field. The sign of the charge 
changes the direction of cyclotron motion so the paths 
are mirrored along Tq. We expect a driving force to break 
the mirror symmetry of the paths. 

In a quantum system this symmetry breaking pertur- 
bation will result in a pseudo-spin dependent phase. We 
can exploit this phase to make an interference measure- 
ment. Specifically, we would prepare the system in an 
initial state |0,f), where |0) is the orbital ground state 
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FIG. 2. The normalized response function for | | for 

the pulse sequence U p , (dashed) or Uc-p, the Carr-Purcell 
like pulse sequence (solid). For both sequences we used At = 
— —r— — = 5? . The frequency of the zeroes between the peaks 
are given by uj = uj±, and the relative peak amplitudes are 
given by e = Note we have scaled the response for 

to the CP pulse sequence by a factor of 16 to account for the 
factor of four increase in interrogation time. 



of the system and <r 2 |t) = It)- Next we place the 
system in a superposition of pseudo-spin states with a 
Ry(n/2) — e -10 "" 71 '/ 4 Raman pulse. Then, we suddenly 
displace the minimum of the harmonic trap by an amount 
ro . If we allow the system to evolve freely in time the two 
different spin states will follow time reversed classical tra- 
jectories. In the process, they also accumulate a pseudo- 
spin-dependent phase term e lcr ^( zxr o) /dtg(t)/i±(t) j wnere 

h±(t) — ^ (uj- s'm(uj + t) — oj + s'm(u)-t)). We now wait 
for a time t at which the two classical trajectories over- 
lap again. We then use a Ry(—n/2) pulse to convert the 
coherence into a population. A S z spin measurement will 
give 



(S,) 



sin 



At [z x r ) • g(i') MO 



(3) 



where we have neglected higher order terms in g/l ui 2 
with the harmonic oscillator length l Q = yjh/mu). We 
can represent the above pulse sequence as the unitary 
matrix U p = R$(— K/2)U(t)D[Rr }Ry(7r/2) where D[r ] 
is a spatial displacement of ro and the time evolution 
operator U(t) can be found exactly. The expectation 
value in Eq. Bis then given by (S z ) = (0,t|Z7tS Jf ?7p|0,t} 
and can be shown to reproduce Eq. [3j 

We now find the response of our interferometer 
to an arbitrary time varying force. We can ex- 
press our spin population measurement as (S z ) — 
sin [/ %g x (u)F (u)] where g(t) = / f^e— 4 g(w) and 
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is the response function of the system, where f(oj) 



sin(^t/2) -iut/2 
jt/2 



The behavior of the response function 
can be seen in Fig. [2] The peak response of the system is 
at the frequencies lu — lu± with relative peak amplitudes 
of lu-/lu + . The bandwidth of the system varies with 1/t 
giving a large bandwidth at small times. Note that our 
system is sensitive to DC signals since F(ui) is finite for 
lu — >• 0. For the purposes of this paper this DC sensitivity 
is unwanted, and we will discuss methods of dealing with 
it later. 

It is important to note that we have waited until the 
coherent states overlap fully. A measurement at a dif- 
ferent time would suppress our signal by a factor of 

A = e~( 2j i h±<,t ^ . The double-exponential suppres- 
sion thus necessitates we obtain as complete an overlap 
as possible. We can eliminate DC signals while simul- 
taneously improving the overlap and eliminating other 
sources of error through application of Rg(n) pulses in 
a sequence analogous to the Carr-Purcell pulse sequence 
in NMR [3T] . If applied at a time when the velocity van- 
ishes, a 7r pulse will time reverse the particle's motion, 
causing it to retrace its path. Such a time is guaranteed 
and occurs at time intervals of t = 27m = ™ . If we 
wait for an additional time time 2t after the first it pulse 
to apply a second n pulse, the path will be time reversed 
again and return to the origin. It is clear that any DC 
signals will be canceled by such a pulse sequence as the 
path returns to itself so the average position is zero. We 
note that this pulse sequence will also help to cancel cer- 
tain noise sources such as Zecman fields or small trapping 
asymmetries. 

In the operator language such a pulse sequence has the 
form 



7T— pulse' 





n— pulses^ 
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FIG. 3. 3(a) The classical path a particle will follow with 



the CP pulse sequence given by Eq. [5] The red path corre- 
sponds to the initial free evolution for a time t. The dashed 
path is the time reversed path allowed to evolve for a time 
2t. Finally, the solid blue path is the return trajectory for 
another t. Note that the three trajectories will overlap in 
practice, and have been offset just for a visual aid. The three 
arcs correspond to the direction of motion of the classical 
trajectory. We have also chosen e = ui+/uj- = 22. 3(b) 
A schematic of the experimental setup. A thermal cloud of 
radius rt is displaced by ro = Ti — r±, which is limited by 
the laser inhomogeneity radius, r;. The two spin states orbit 
along trajectories mirrored across ro. 3(c) The dependence of 



the sensitivity of the system based upon atom number. Be- 
low N c ~ 10 6 the sensitivity grows as . 1 Above N c the 

sensitivity S ~ 10 _7 m/s 2 /v / Hz is independent of the number 



U C p = Ry{-n/2)U(t)Ry(ir)U(2t)Ry(Tr)U(t)D[Rr ]Ry(Tr/2). oi particles. |3(d)| The bandwidth of the system with optimal 

(5) sensitivity as a function of the number of particles. 

Application of such a pulse sequence modifies the re- 
sponse function to 



F(u) = 2isin(wi) [F (u))e iut + F£(u)e 



■iu)t~\ Alut 



(6) 



where Fo(u>) is the response function given in Eq.|4| The 
complex conjugate term Fq (lu) arises due to the time 
reversal of the paths. This response function is plotted 
in Fig. [2j The new response function now vanishes at 
lu = and lu = lu± , however we still have large sensitivity 
in the frequency range u = | around lu± . 

We now generalize from a single particle to a ther- 
mal ensemble of ultracold dilute atoms in the presence 
of an induced spin-orbit Hamiltonian of the form in 
Eq. [T] From a practical standpoint a system of di- 
lute cold atoms allows us to neglect complications aris- 
ing atom-atom interactions in a BEC. Consider a ther- 
mal ensemble of cold atoms at a temperature T. Using 
the Glauber P-representation [32] the density matrix has 
the form p = J da e -l Q +l 2 /<™+)-|a-l 2 /<™-> | a ) ( a | ; wne re 



(m) = [exjp[hu!i / kT] — is the average occupation 
for the classical mode of frequency w^. Such an ensem- 
ble suppresses the expectation value of S z operator by 
a factor e _ ^"+^ 7 +' < rl — > I t — I relative to the single parti- 
cle/zero temperature expectation value. The suppression 
factor 7± depends on the pulse sequence used. For ex- 
ample, using the pulse sequence U p given above we get 
7+ = l f dt< (g x + igyY^' and 7 _ = l f f< dt< (g y + 
*ffx)e IW_ ' . Note that we can express J dt r • g as a super- 
position of 7+ and 7_. For the Carr-Purcell like pulse 
sequence the suppression factor j± is more complicated, 
but can still forms a basis for which we can express the 
phase. This implies that j± has a similar frequency de- 
pendence to F(lu). 

We are now in a position to discuss the measurement 
capabilities of such a system. We first estimate the max- 
imum AC signal such a system can measure. To avoid 
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signal suppression due to the finite temperature of the 
ensemble we bound the maximum strength of the AC 
signals our system can measure by g max < "J^y mni " ^ n 

this limit the sensitivity for our detector can be estimated 
with 

where the lifetime of one measurement, 1/r = 7 se + 7 co n 
is limited by spontaneous emission, j se and collisions, 

^coll • 

To minimize the sensitivity of the system we need to 
consider the effect of collisions and the spatial configu- 
ration of the system. We desire to confine our system 
to a "laser homogeneity" radius, r;, for which nonlinear 
variations in the laser fields are suppressed. The effec- 
tive 2D system will use an axial trapping potential of 
w || > \og(2)kT/h to freeze all motion into a single trans- 
verse mode. Thus our system will have Ni = ri/d layers, 
where d = yjhjmuj^ , providing for an increase in sen- 
sitivit of \j\jN\. The radius n further constrains our 
sensitivity by bounding the maximum trap displacement 
by r o — r i ~ r t, where r t — \J {n)l ~ (v) jib is the ther- 
mal radius of a thermal ensemble and (n) — kT/hui and 
(v ) = yJikT jm are the respective high temperature ther- 
mal occupation number and velocity. (Fig. |3(b)| 

The lifetime of the system will be dominated by spon- 
taneous emission at low densities and collisions at high 
densities. To optimize the sensitivity we desire to place as 
many atoms per layer as possible. The collisional scat- 
tering rate is given by 7 co ;/ = Na ^l a , where a is the 
interparticle scattering length. The critical number of 
atoms at which the collision rate begins to dominate the 
spontaneous emission rate is N c — ■ We see that in 

the small atom number limit the sensitivity is a monoton- 
ically decreasing function of the trapping frequency and 
has a l/\/N a dependence. However, in the large atom 
limit the sensitivity is minimized at a trapping frequency 
of Uimin = 2(v)/ri and the sensitivity becomes indepen- 
dent of the number of atoms per layer. Note that in this 
limit the bandwidth of the system is increased by adding 
atoms. (See Fig. [3(d)] ) 

We assume our thermal ensemble has a temperature of 
T ~ lmK with a frequency scale a) = 2-/rkHz. At these 
temperatures the gas is non-degenerate and is described 
well by a classical gas. For this temperature we find an 
upper bound of g ~ 10~ 2 m/s 2 before exponential sup- 
pression of the signal above becomes relevant. We will 
consider a cold gas of 87 Rb cooled to T — 1/iK with an 
axial confinement distance of d — 1/im. We take the 
spontaneous emission rate to be r se = l/70ms |23| and 
laser inhomogeneity radius to be ri = 10 — 25/im. In 
the N a ^> N c limit we estimate the sensitivity to be 
S ~ 10 -7 ^^^. A similar analysis for a 3D system gives a 



sensitivity drop of approximately an order of magnitude. 
We note that had we instead used a fermionic species 
we would obtain a similar result since we have two spin 
species. 

The concept of a continuous coupling of spin to mo- 
mentum can also be extended to a continuous coupling 
of spin and position. We note that in a harmonic trap 
position and momentum are dual variables, and thus a 
spin-dependent term in the Hamiltonian that has spatial 
variation will experience a similar phase accumulation 
to the system described above. An example would be 
trapped spin-1 system in the presence of a Zeeman field 
with a strong spatial variation. In such a system the 
Zeeman field will act to trap (anti-trap) the S z = +1 
and S z = — 1 spin states with different trapping poten- 
tials. This will play a similar role to the opposite charge 
couplings to gauge fields given above. However, such a 
system requires strong magnetic field gradients so it may 
be impractical. 

Finally, we note that this system is not limited to mea- 
surements of AC signals. Through appropriate modifi- 
cations of pulse sequences such a scheme is capable of 
measurements of both gravitational gradients as well as 
constant DC gravity. For example, the system will be 
sensitive to DC signals if a measurement is made for any 
path at which the average position perpendicular to the 
axis of displacement is non-zero. Since overlap is guar- 
anteed at the anti-node of the orbit, this would be a 
reasonable measurement point. Similarly, the system is 
sensitive to gravitational gradients if a measurement is 
made at a point of overlap when no ir pulses have been 
applied. Both of these effects will be eliminated through 
a Carr-Purcell like pulse sequence. Due to electronics 
noise, the sensitivity of our system to either DC or gra- 
diometric signals will be significantly lower than existing 
setups. Thus, we have chosen to not pursue them in de- 
tail. 
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